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Abstract 



I discuss three different selected topics in top quarlc physics. The first topic concerns the 
next-to-next-to-leading order calculation of the hadroproduction of top quark pairs and the 
role of multiple polylogarithms in this calculation. I report on an ongoing next-to-next- 
to-leading order calculation of heavy quark pair production in hadron collisions where the 
loop-by-loop part of the calculation is about to be completed. Calculating the loop-by- 
loop part allows one to take a glimpse at the mathematical structure of the full NNLO 
calculation. The loop-by-loop contributions bring in a new class of functions introduced 
only eight years ago by the Russian mathematician Goncharov called multiple polyloga- 
rithms. The second topic concerns a next-to-leading order calculation of unpolarized top 
quark decays which are analyzed in cascade fashion t — > b + W~^ followed by —>■ + 
Finally, I present some next-to-leading order results on polarized top quark decays which 
are analyzed in the top quark rest system. 



1 Introduction 



1 begin my talk with a few remarks on present top quark yields at the Tevatron and 
on expected top quark yields at the LHC which will start running at the end of 2007. 
After a slow start in early 2001 Tevatron II started reaching design peak luminosities of 
8.5 • 10'^^cm~^s~^ in 2004. The best weekly performance was early 2006 with a weekly 
integrated luminosity of 25pb~^. If Tevatron II could perform at this rate it would be able 
to collect 1.3 fb~^ in a year. There has been a three months shutdown in the spring of 2006 
with some (electron cooling) improvements on the p beam. The hope was that there will be 
a factor two or three improvement in luminosity after the shutdown. Such a factor would 
be dearly needed if one wants to reach the projected total of 8fb~^ when the machine is 
closed down in 2009. At the time of writing this factor has not been realized so far after a 
few months of post-shutdown running although the machine is performing quite well with 
continuous improvements. At a cm. energy of ^/s = 1.96 TeV with a{tt) ^ 6.8 pb one 
expects around 7000 ti pairs at each detector (CDF and DO) for an integrated luminosity 
of 1 fb~^. Single top production occurs at about 33% of the ti pair production rate but has 
not been detected so far. 

Much bigger samples of top quarks will be available at the LHC. Due to the higher 
energy of the LHC the cross section increases by a factor of 100. Also there will be a ten- 
fold increase in luminosity at the LHC. Thus one will have 10'' tt-pairs per year, or one 
tt-pair every four seconds at each detector (ATLAS and CMS). In a later high luminosity 
run there will be another factor of ten increase in luminosity such that one will have a 
tt-pair produced every half second. Again single top production occurs at approximately 
one-third the rate of tt-production. Singly produced top quarks will be highly polarized 
because they are produced weakly. This opens the way to study angular correlations 
between the polarization of the top quark and its decay products which forms the third 
topic of this talk. 

The yield of top quark pairs at the International Linear Collider (possibly starting in 
2015) will be ^ (1 — 4) ■ lO^/y depending on the cm. energy ?a 360 — 800 GeV. Ine"*" — e~- 
interactions a high degree of polarization of the top (or antitop) quark can be achieved 
through tuning of the beam polarization. 

2 NNLO description of heavy top quark production 

The full next-to-leading order (NLO) QCD corrections to hadropro duct ion of heavy flavors 
were completed as early as 1988 [HE]- They have raised the leading order (LO) estimates 
[3j but were still below the experimental results on bottom quark pair production (see 
e.g. 0). In a recent analysis theory moved closer to experiment j4i- First experimental 
results on hadronic tt-pair production Ej are in agreement with theoretical NLO QCD 
predictions [3 |H] within the large theoretical and experimental error bars. 

A large uncertainty in the NLO calculation results from the freedom in the choice 
of the renormalization and factorization scales. These scale uncertainties amount to a 
~ 10% theoretical error in the NLO cross section predictions [7]. The dependence on the 
factorization and renormalization scales is expected to be greatly reduced at next-to-next- 
to- leading order (NNLO). This will reduce the theoretical uncertainty. 
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a) 




d) 

Figure 1: Exemplary gluon fusion diagrams for the NNLO calculation of heavy hadron 
production 

In Fig. n we show one generic diagram each for the four classes of gluon induced contri- 
butions that need to be calculated for the NNLO corrections to hadropro duct ion of heavy 
flavors. They involve the two- loop contribution (Fig. the loop- by-loop contribution 
(Fig. IHd), the one- loop gluon emission contribution (Fig. [TJ:;) and, finally, the two gluon 
emission contribution (Fig. ^). In our work we have concentrated on the loop-by- loop 
contributions exemplified by Fig. ^p. Specifically, working in the framework of dimen- 
sional regularization, we have calculated O(e^) results for all scalar massive one-loop one-, 
two-, three- and four-point integrals that are needed in the calculation of hadronic heavy 
flavour production JU]. Because the one- loop integrals exhibit infrared (IR)/collinear (M) 
singularities up to 0{e~^) [H] one needs to know the one-loop integrals up to C(e^) since 
the one-loop contributions appear in product form in the loop-by-loop contributions. It is 
exactly the O(e^) terms in the scalar massive three- and four-point integrals that bring in 
multiple polylogarithms |TU1 ITT]. 

Calculating the loop-by-loop contributions allows one to obtain a glimpse of the com- 
plexity that is waiting for us in the full NNLO calculation. This complexity does in fact 
reveal itself in terms of a very rich polylogarithmic structure of the Laurent series expan- 
sion of the scalar one-loop integrals as well as the appearance of multiple polylogarithms 
of maximal weight and depth four. 

To underscore the statement that the loop-by-loop contributions reveal the mathemat- 
ical structure of a full NNLO calculation let us take a look at the paper by Bernreuther 
et al. ^21 who calculated the 0{e^) contributions to the one-loop vertex correction of the 
process V —>■ QQ. These are needed for the loop-by-loop part of a NNLO calculation of 
heavy quark pair production in e"*" — e~-annihilations. The result can be expressed in terms 
of one-dimensional harmonic polylogarithms of maximal weight three. The same paper 
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also lists results on the corresponding 0{e^) two-loop vertex corrections which contain one- 
dimensional harmonic polylogarithms of maximal weight four. This shows that the same 
mathematical complexity appears in the loop-by-loop contributions as in the two-loop 
contribution. Let me mention that heavy quark pair production in — e~-annihilations 
is a somewhat simpler problem than heavy quark pair production in hadronic collisions 
because of the appearance of one additional mass scale in the latter case. This explains 
why one has only one-dimensional harmonic polylogarithms in — e~ -annihilations case 
compared to the multiple polylogarithms appearing in the hadronic collision calculation. 
Even then the two-loop vertex correction to V ^ QQ listed in ^2] takes up more than 
twelve pages. 

The scalar four-point integrals appearing in the calculation of the loop-by-loop eval- 
uation are the most difficult to calculate. They contain a very rich structure in terms of 
polylogarithmic functions. For example, the e^-coefficients of the Laurent series expansion 
of the four-point integrals contain logarithms and classical polylogarithms up to order four 
{i.e. Li^) in conjunction with the C^functions C(2),C(3) and ^(4) and products thereof, 
and a new class of functions which are now termed multiple polylogarithms [13 . 

Since this is a conference on mathematical physics it is appropiate to dwell a little on 
the subject of multiple polylogarithms. A multiple polylogarithm is represented by 



The indices and k are positive integers. The multiple polylogarithms are classified 
according to their weight w = mi + m2 + ... + and their depth k. We mention that a 
very efficient program for the numerical evaluation of multiple polylogarithms has recently 
been developed in Mainz which, characteristically, is based on the language GiNaC [1^. 

The classical polylogarithms, Nielsen's generalized polylogarithms, the one- and two- 
dimensional harmonic polylogarithms are all special cases of Goncharov's multiple poly- 
logarithms (see e.g. fH])- For example, the classical polylogarithms 



Li, 




where the iterated integrals are defined by 




o ... o 



dt 





(1) 



are multiple polylogarithms of weight n and depth 1. 
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In our original Feynman parameter calculation our results were written down as one- 
dimensional integral representations given by the integrals 



p . N /"^ J + ^^y) l^("2 + cr2l/) ln(a3 + CT^y) 

-^'<7i,72a3(ai,a2,a3,tt4) = / dy ■ (2) 

^0 cxA + y 



and 

"i ln(ai + criy)Li2(a2 + "si/) 



F^j(ai,a2,a3,«4) = / "Z/ ^ , (3) 

JO 0!4 + y 

where the cxj take values ±1 and the a/s are combinations of the kinematical variables of 
the process. The numerical evaluation of these one-dimensional integral representations 
are quite stable. The functions F^.^ 0-20-3 «2, «3, ^4) and F^^{ai, 02, 03, 04) are related to 
multiple polylogarithms of maximal weight and depth four as shown in [TT] . 

We are now in the process of computing the full loop-by-loop contributions including 
the spin and colour algebra arising from squaring the full one-loop amplitudes as given in 

. A first result has been obtained for the Abelian case of photon-photon production of 
heavy quark pairs [T7] . 



3 Decays of unpolarized and polarized top quarks 

After this brief mathematical detour I return to the physics of top quark decays. In the SM 
the top quark decays almost 100% to a and a bottom quark. Also, the top quark decays 
so fast that it retains its initial polarization when it decays. I describe both unpolarized 
and polarized top quark decays. In the unpolarized case I analyze top quark decays in 
cascade fashion as a two step process involving the decay t b + and + ui 

in the respective rest frames of the top quark and the VT^-boson. In the polarized case 
I perform the decay analysis in the rest system of the top quark. I discuss polar and 
azimuthal correlations involving the polarization of the top quark and the momenta of the 
decay products in the decay t(t) ^ Xf, + + i//. 

The decay of a polarized top quark into a H^+-boson and a jet with 6-quantum numbers 
t(t) — ^ Xft + /+ + z/; is desribed by altogether eight invariant structure functions (see e.g. 

HHiiini). 

- {q-st){- g''" Gi + pf^r G2 - le'-'^'-pt^.q^ Ga) + (4) 
+ (sr^r + s^V) Ge + le^'^^'ptpsta Gs + ie'''""'q,sta Gg , 

There are three unpolarized structure functions Hi 2,3 and five polarized structure functions 
from the set Gi^2,3,6,8,9 ^- In general the invariant structure functions are functions of 
go and q^. In the narrow resonance approximation for the Vr"'"-boson, which we shall 
adopt in this talk, one has = m^. The aim of the game is to measure the different 
unpolarized and polarized structure functions (or moments thereof) and to compare them 

^In physical expressions the three structure functions G3,Gs and Gg contribute only in two pairs of 
linear combinations 
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to theoretical predictions. The different structure functions can be separated since they 
contribute to the rate with different dependencies on the electron energy and, in the case 
of the polarized structure functions, they can be measured through polar and azimuthal 
correlations involving the polarization direction of the polarized top quark. 



Unpolarized top quark decays 

In the decay t ^ Xj, + the is polarized. The is self-analyzing in the sense 
that the angular decay distribution of its decay products — >• l'^ vi can be used to 
reconstruct the polarization of the W^. We shall analyze the unpolarized decay in cascade- 
type fashion, i.e. we shall analyze the decay —>■ -\- ui in the rest frame of the 
W^. This brings in the three unpolarized helicity structure functions Hjp^, Ht^, Hl (or for 
short H^, H_, Hl) which are linearly related to the three unpolarized invariant structure 
functions Hi,H2, via ^Hl 

H+ = Hi + \q\mtH3, (5) 
if_ = Hi-\q\mtH3, (6) 
Hl = ml,H^ + \q\^mlH2, (7) 

The polar angle decay distribution is given by 

i/^- = |(1 + cos^)^ + |(1 - cosef H^ + \ sin' 9 . (8) 
r a COS 6* 8 8 4 

where the angle 6 is defined in Fig|21 The 7^+,7Y_ and Tii are the normalized transverse- 
plus, transverse-minus and longitudinal helicity structure functions, resp., such that + 

n- + nL = i. 

From the polar angle dependence in Eq.(jH} or from matching m-quantum numbers in 
the rest frame decay (see Fig|21) it is clear that 

7^+ : favours forward 
Ti- : favours backward 

Translated to the top quark rest frame this implies that JF+ (JF_) produce harder (softer) 
Z"^'s which can be used to experimentally separate the contributions of the three helicity 
structure functions. 

At the Born term level the SM prediction is (m^ = 175 GeV, = 0) 



7^+(Born) = (forbidden) (9) 

n-{Born) = = 0.297 

H^Born) = = 0.703, 

where y = mw/'fnt- At the Born term level, with mi, = 0, 7i+ is not populated because 
of angular momentum conservation in the two-body decay process t b + where for 
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Figure 2: Definition of tlie polar angle 9 in tlie rest frame decay of + ui. The two 

lines "/ /" indicate a boost to the rest system of the . The arrows next to the lepton 
lines give the helicities of the leptons. 



mft = the bottom quark has 100% negative helicity and the 6-quark and the are in 
a back-to-back configuration. 

The present experimental results on Ti^ are consistent with zero within large error bars. 
For example, using 230 pb^^ DO quotes a value of ?-^+ = 0.00 ± 0.13 (stat) ± 0.07 (syst) 
[20] • CDF finds n+ = 0.00 +[|;3J (stat + syst) or n+ < 0.27 at the 95% confidence level 
fIJ\ . Using the same data sample of 200 pb^^ CDF quotes a value of Hl = ^■'^^^034 fo^ 
the longitudinal helicity of the H^"'"-boson, also compatible with the SM prediction. 

The vanishing of 7i+ is no longer true for additional gluon or photon emission, or when 
one takes into account bottom mass effects. When all of these are taken into account one 
has 1221 ES! 

n+ = 0.00102(QCD) + 0.00008(EW) + 0.00039(m6 ^ 0), (10) 

where the numbers give the 0{as) QCD corrections, the 0{a) electroweak corrections 
and ruh 7^ corrections (mf, = 4.8 GeV). Numerically the correction to ?i+ occurs only at 
the pro mille level. It is safe to say that, if top quark decays reveal a violation of the SM 
(V — A) current structure that exceeds the 1% level, the violations must have a non-SM 
origin. 

The results for the corresponding corrections to Ti- and Tii are listed in terms of rates 
normalized to the total Born term rate, i.e. F, = ri/r(Born). The normalized partial 
Born term rates rj(Born) are factored out. Corrections coming from NLO QCD, from the 
NLO electro- weak corrections (EW), from the finite width correction (BW) and from 
rrih 7^ effects are listed separately. One has 

f_ = 0.297 [1 - 0.0656(QCD) + 0.0206(EW) -0.0197(BW) - 0.00172(m;, ^ 0)] 
Tl = 0.703 [1 - 0.0951(QCD) + 0.0132(EW) -0.0138(BW) - 0.00357(mfe ^ 0)] 

Written in terms of the normalized Tii this translates into a +2.4% upward shift from 
7f_(Born) = 0.297 and a -1.2% downward shift from ^^(Born) = 0.703. Judging from 
the fact that Tii and 7i_ will eventually be measured to better than 1% it is quite clear 
that one has to take radiative corrections into account when comparing experiment with 
theory. 
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Figure 3: Top mass dependences of the ratio Tii = r^^/F. Full line : Born term. Dashed 
line: Corrections including (QCD), electroweak (EW), finite- width (FW) and {nib ^ 0) 
Born term corrections. 



In Figiniwe show the top mass dependence of the ratio Tii = Ti/r- The horizontal 
displacement of the Born term curve and the corrected curve is ~ 3.5 GeV. One would thus 
make the corresponding mistake in a top mass determination from the measurement of Tii 
if the Born term curve was used instead of the corrected curve. If one takes uit = 175 GeV 
as central value a 1% relative error on Tii would allow one to determine the top quark 
mass with an error of ~ 3 GeV. 



4 Polarized top quark decays 

Contrary to the analysis of unpolarized top quark decays described in the last subsection 
polarized top quark decay will be analyzed altogether in the rest frame of the decaying 
top quark. This is the natural choice for an experimental analysis. Choosing a particular 
two-particle rest subsystem is only of advantage if that particular subsystem is resonance 
dominated as was discussed in the unpolarized decay case. 

The general angular decay distribution of the rest frame decay of a polarized top quark 
decaying into a jet and a lepton /+ and a neutrino is given by 



ax idqod cos opa(p Air \axidqo axidqo dxidqo J 

where the polar and azimuthal angles Op and describe the orientation of the polarization 
of the top quark relative to the decay plane formed by the decay products of the top 
quark. The scaled energy and the scaled mass of the are denoted by go = qo/'^t and 
y = mw/mt. As usual we define a scaled lepton energy through xi = 2Ei/mt. P is the 
magnitude of the top quark polarization. F^ stands for the unpolarized rate, and F^ and 
Tc stand for the polar and azimuthal correlation rates. In [22] we have considered three 
different helicity systems to analyse the polar and azimuthal correlations in the rest frame 
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(la) 




(3a) 

Figure 4: The definition of tfie polar angle Op and tlie azimutlial angle in tlie rest frame 
decay of a polarized top quark in three different helicity systems. The event plane defines 
the (x, z)-plane with (la) pi \\ z and {py)x > 0, (2'a) q \ z and {j)])x > 0, and (3a) || z 
and {:pi)x < 0. 



decay of a polarized top quark as shown in FigEI 

It is important to realize that correlation measurements in each of the helicity frames 
constitute independent measurements of the invariant polarized structure functions. To 
illustrate this point let us consider the contribution of the invariant polarized structure 
function Gi to the polar and azimuthal correlations in the above three helicity systems. 
The decay rate is proportional to L^^H^j^. One then obtains 



mo - y 



ql - ?/2 cos Op2 

qo-y^ n 

COSfc/ps - 



I 2qo - xi 



y 



cosOpi H \/xi(2qo — xi) — y"^ sin 6'pi cos( 

Xi " 



y 



2qo - xi 



Xi{2qo -xi) -y"^ smOp^ cos. 



(12) 
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where the three contributions in the curly bracket refer to the polar and azimuthal corre- 
lations in the three helicity coordinate systems with the 2;-axes along (1) the lepton , 
(2) the VT^-boson and (3) the neutrino ui . From Eq.()12|) it is clear that Gi contributes 
quite differently to the correlation functions in the three reference systems. 

In j2n] we have calculated the Born term and NLO QCD contributions to the polar 
and azimuthal correlation functions dTs and dTc in the three different helicity systems. 
We were able to obtain closed form expressions for the totally integrated angular decay 
distributions. The results are too long to be listed here but can be found in [23]. We 
mention that we find agreement with [23 EH IIHI IHOj for the unpolarized case (ITa and the 
polar correlation function dTs in systems 1 and 3. In numerical form one has 
z-axis along lepton (system (la)) 



^-pNLO p{0) 



A 



dcos6pd(j) Att l 
2;-axis along IV """-boson (system (2'a)) 



(1 - 8.54%) + (1 - 8.72%)Pcos0p - 0.24% P sin cos (13) 



d COS 9pd(f) All 
^-axis along neutrino (system (3a)) 



1 - 8.54%) + (0.406 - 11.62%)P cos dp - (0.760 - 8.20%) P sin Qp cos 

(14) 



d cos 9pd(f) 4:71 



[1 - 8.54%) - (0.318 - 1.02%)Pcos^p - (0.919 - 8.61%) P sin ^p cos 

(15) 

In all the three expressions we have factored out the Born term rate F^"* . The first number 
in the round brackets stands for the LO Born term rate whereas the second number gives 
the percentage change due to the NLO QCD corrections. 

Let me first discuss the LO correlation functions. I shall refer to F^/F^ and F^/F^ 
as the polar and azimuthal analyzing power, respectively. In system (la) along z) 
the polar analyzing power is 100% which necessarily implies that the azimuthal analyzing 
power is zero in this system. In fact, the vanishing of Fc in system (la) can be seen to 
directly follow from the left-chiral {V — A) structure of the SM quark and lepton currents 
[2n]- The polar analyzing power in the systems (2'a) and (3a) is less than 100% with +41% 
and -32%, respectively. As mentioned before the LO azimuthal analyzing power in system 
(la) is zero. In system (2'a) and (3a) the azimuthal analyzing power is reasonably large 
with -76% and -92%, respectively. 

Except for the polar correlation in system (3a) all NLO corrections go in the same 
direction. They reduce the LO results by approximately 10%. This implies that the polar 
and azimuthal analyzing powers are not changed very much from their Born term values 
through radiative correction. An exception is system (3a) where the polar analyzing power 
is changed from -31.8% to -34.4%. This amounts to a 8.2% change in analyzing power 
through radiative corrections which is surprisingly large. 
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5 Summary and conclusions 

In this talk I have covered three selected topics in top quark physics. The first topic con- 
cerned the NNLO calculation of hadronic top quark pair production where the loop-by- 
loop part is now being completed. The other three missing parts of the NNLO calculation 
(two-loop, one-loop gluon emission, two-gluon emission) are more difficult and will very 
likely take another five to ten years to complete. Such a large-size calculation will require 
a dedicated international effort of the theoretical community which will have to be coordi- 
nated by one of the big international centers of particle physics. In the second and third 
topic I discussed NLO QCD predictions for unpolarized and polarized top quark decays 
which should be amenable to experimentals tests in the next coming few years. 

Acknowledgements: I would like to thank my collaborators S. Groote, W. S. Huo, 
A. Kadeer, D. Kubistin, Z. Merebashvili and M. Rogal for participating in the work I 
have reported on in this talk. My thanks are also due to Riazuddin and F. Hussain for 
organizing such a wonderful conference in Islamabad, and to all the Pakistani graduate 
student helpers who made the meeting such a joyful affair. 
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